function [] = system1(system_no,r,L,M,a,b)
    fprintf('** System no %d\n',system_no);
    disp('* Most positive representable number');
    largest_representable_number(r,L,M,b);
    disp(' ');

    disp('* Most negative representable number');
    largest_neg_representable_number(r,L,M,a);
    disp(' ');

    disp('* Number of representable numbers');
    nrn = ((L+M)*(b-a)/2+1)^(L+M);
    fprintf('((L+M)*b+1)^(L+M) = ((%d+%d)*%d+1)^(%d+%d) = %d\n',L,M,b,L,M,nrn);
    disp(' ');

    disp('* ulp');
    fprintf('r^-L = %d^-%d = %.32f\n',r,L, r^-L);
    disp(' ');
    
    disp('* Redundancy index');
    p = -a+b-r+1;
    fprintf('p = -a+b-r+1 = -%d+%d-%d+1 = %d\n',a,b,r,p);
    disp(' ');
end

